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Abstract 

Super symmetric SU(5) GUT augmented with anomaly free U(1)f flavor symmetry is 
presented. Very economical field content and U (1) f charge assignment are obtained by specific 
construction. In particular, three families of 10+5 chiral matter, along the SU (5) singlet states 
(some of which serve as right handed neutrinos) are obtained. Appealing texture zero Yukawa 
matrices provide natural understanding of hierarchies between charged fermion masses and 
mixings. The model predicts inverted hierarchical neutrino mass scenario with interesting 
implications. 



■ 1 Introduction 

_cp_- Although being very successful, the Standard Model is unable to resolve some puzzles. Among 
them is a problem of fermion flavor. The origin of hierarchies between charged fermion masses and 
CKM mixing angles is unexplained. Is there any underlying theory which might generate these 
hierarchies in a natural way? Moreover, in order to explain the neutrino data [1-3] some extension 
of the Standard Model, generating neutrino masses and mixings, is necessary. The number of 
fermion generations is a mystery. Do we have only three chiral families of quarks and leptons? Is 
any selection rule dictating the number of fermion generations? 

Motivated by these questions, in this paper we address these issues within the framework of su- 
persymmetric (SUSY) SU(5) Grand Unified Theory (GUT). Latter's motivation is to have unified 
description of electro- weak and strong interactions [4], while SUSY provides natural understanding 
of gauge hierarchy problem as well as the successful gauge coupling unification [5] . For understand- 
ing hierarchies between fermion masses and mixings, we apply Abelian flavor symmetry [6] U(l)p 
which by requirement is non-anomalous. The U(l)p, in combination with SUSY SU(5) GUT, due 
to anomaly constraint allows only three chiral families of matter (10 + 5)-plets and few SU(5) 
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singlet states. We use some of these singlet states as right handed neutrinos, in order to build re- 
alist neutrino sector. Anomaly constraints fix U(l)p charge assignment in such a way that texture 
zero quark and lepton Yukawa matrices are generated. Together with natural understanding of 
hierarchies, model predicts inverted neutrino mass hierarchical scenario blending well with recent 
Daya Bay observation [2]. 

2 Three Family SUSY SU{b) x U(1) F 

Consider SUSY SU(5) GUT augmented with anomaly free U(l)p flavor symmetry. Setup with 
anomaly free U(l)p, will alow to gauge U(l)p and remain within conventional 4-dimensional field 
theoretical framework (without need of discussing U (l)'s of a stringy origin [7]). In a recent work [8] 
the way of finding anomaly free U(1)f flavor symmetry within SUSY S77(5) GUT was suggested. 
The finding was realized by embedding of SU(5) x U(1)f in a single non-Abelian group with 
anomaly free field content. 2 In this way, the U(l)p charge assignment can be fixed. Amongst 
several assignments, found in [8], there is one which also dictates the number of generations to be 
three. This, As will be shown below, leads to very economical and attractive scenario for fermion 
masses and mixings. Before showing this, we briefly discuss the way of finding of such U(l)p- 

The states non-trivial under SU(5) group, we introduce, will be just those of minimal SUSY 
S77(5). These are scalar superfields 11(24), 77(5), 77(5) and three families of matter (10 + 5) su- 
permultiplets. We assume that £ is not charged under U(l)p and thus does not contribute to 
the anomalies. Therefore, upon finding anomaly free U(1)f charge assignment we will deal with 
three 10-plets, one 5-plet (which is 77) and four 5-plets (=three matter 5-plets plus 77). As already 
mentioned, we search U(1)f charge assignment by embedding of SU(5) x U(1)f in non-Abelian 
group with anomaly free matter [8]. Let us consider SU(7) group with chiral supermultiplets 
35 + 2 x 7. This simple set is anomaly free [10]. Here 35 is three index antisymmetric repre- 
sentation and 7 is an anti-fundamental of SU(7). Decomposition of these states via the chain 
SU(7) ->■ SU(6) x U(l) 7 ->■ SU(5) x U(l) 6 x U(l) 7 looks 

35 = 20 3 + 15_ 4 = (10-3 + 10 3 )3 + (10 2 + 5_ 4 )_ 4 , 

7 = 6_! + 1 6 = (5_! + l 5 )-i + (lo)e , (1) 

where subscripts inside and outside of parenthesis indicate U(1)q and U(l)j charges respectively. 
Note that U(1)q and U(l)j are coming from SU(6) and SU(7) respectively, with corresponding 
generators Y u(lh = ^Diag (1, 1, 1, 1, 1, -5) and Y u(1)7 = ^fDiag (1, 1, 1, 1, 1, 1, -6). The nor- 
malization factors ^= and ^= are omitted in Eq. (1). Now, in (1), without change of charge 
assignments, we replace the pair of (10 + 5)-plets by the pair (10 + 5). With this replacement all 
anomalies (at the level of S77(5) and C/(l)'s) will remain intact, i.e. will still vanish. Thus, we will 
have the following anomaly free content 

(10_ 3 + 10 3 ) 3 + (10 2 + 5_ 4 )- 4 + 2 x [(5_i + l 5 )-i + (l )e] , (2) 

which involves three families (!) of matter (10 + 5) supermultiplets plus four 577(5) singlets. Some 
of these singlets will be applied as right handed neutrinos (RHN). Worth noting that, in difference 

2 The way of this finding differs from those used earlier [9]. 
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from 5*0(10) GUT, the £77(5) does not involve (require) the RHN states. In the £77(5) scenario, 
we have just built up, the RHNs are required for anomaly cancellation. 

By Abelian symmetries U(l) e and U(l) 7 , with charges given in Eq. (2), we can build superpo- 
sition aQu(i) 6 + bQu(i) 7 - This superposition is automatically anomaly free for arbitrary a and 6, 
because the orthogonal generators 5 / f/(i) 6 an d ^7/(i) 7 originate from single 577(7). Thus, using (2) 
we can write the anomaly free set 

10-3a+36 + 103a+3b + 102a-4ft + ^>-Aa-4b + 2 X (5_^-6 + I5S-6 + ^6b) ' (^) 

where subscripts denote charges. As it turns out, for building realistic phenomenology it is useful 
to add to this superposition another anomaly free U(l), which can be found by similar procedure. 
For instance, consider 27-plet of E 6 group with a chain E 6 -» £0(10) x U(1) E6 ->■ SU(5) x C/(l) Eg 
of decomposition [11]: 

27 = 16x + 10_ 2 + 1' 4 = (10 + 5 + l)i + (5 + 5')_ 2 + 1 4 • (4) 

Here subscripts denote £/(1)b 6 charges. In this content, we can replace 5_2 with 5_2 and at the 
same time add 10 p + 10_ p . Moreover, we can add two SU(5) singlets with U(1)e 6 charges k and —k 
respectively. With this replacement and additions, the anomalies ££7(5) 3 , SU(5) 2 -U(1)e 6 , etc, will 
be unchanged. This content allows to build superposition of three Abelian groups U(l)e, U(l) 7 and 
U(1)e 6 - Qsup = a<5f/(i) 6 + bQu(i) 7 + cQu(i) E6 . Thus, the field content and Q sup charge assignment 
will look: 

10-35+35+pc + ^3a+3b~pc + ^2a-4b+c + ^-Aa-Ab-2c + ^-a-b+c + ^-a-b-2c 

+h- a -l + - c + l 5s -- 6+ 4 C - + l' 66+ fcc + i'eS-fce > with 30 «( 3 + 2 P) = c(2A; 2 + 10p 2 - 27) . (5) 

Relations between a, c, fc and p (imposed for c ^ 0) given in Eq. (5) insures that all anomalies 
vanish. Clearly, with rational selection of a, k and p the value of c also will be rational. The set 
given in Eq. (5) is one simple selection among several options and opens up many possibilities for 
model building with realistic phenomenology. We will identify these charges with the charges of 
U(1)f flavor symmetry. 

3 Model: Quark and Charged Lepton Yukawa Textures 

To the field content of Eq. (5) we add the pair 5 q + 5_ g . This is needed to have, besides the matter 
fields, the Higgs supermultiplets H + H. Thus, in total we have three 10-plets, one 5-plet and 
four 5-plets plus SU(5) singlets. The U(1)f charge assignment (q, —q for 5 q , 5_ g and for remaining 
5-plets given in Eq. (5)) is not unique. We can exchange 5-plet 's U(1)f charge with one of the 
5-plets' charge. With this, all anomalies will still vanish. In addition, out of the four 5-plets, any 
of them can be identified with the Higgs superfield H. As it turns out, for the charges of the pair 
(H,H), we will have 13 possible options for (Qh,Qh) : 

(QH,Qs) (l) ={(q,-q), (g, -4a -46 -2c), (q,-a-b + c), (q,-a-b-2c), 

(_4a - 46 - 2c, -a - b + c), (-4a - 46 - 2c, -a - b - 2c), (-4a - 46 - 2c, q), 

(—a — 6 + c, —4a — 46 — 2c), (—a — 6 + c, —a — 6 — 2c), (—a — 6 + c, q), 

(-a - 6 - 2c, -4a - 46 - 2c), (-a - 6 - 2c, -a - 6 + c), (-a - 6 - 2c, q)} , (6) 
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with I — 1, • • • , 13. Note that we have left out possibilities obtained from those given in (6) by the 
substitution q — > —q. These 13 options open up various possibilities for the model building [12]. 
Below we present one of them, which we found to have nice and attractive properties with interesting 
implications for fermion masses and mixings. 

U(1)f symmetry breaking 

In order to break U(1)f gauge symmetry, we introduce SU(5) singlet pair of flavon superfields 
X + X with U(1)f charges 

Q(X) = -1 , Q(X) = 1 . (7) 

Without loss of generality, we have normalized Savons' charges modulo to one. The scalar compo- 
nents of X and X acquire VEVs 3 

(\X\)_ (\X\)_ 

- e > ~n~ - 6 > I 8 ) 



M P1 Mpi 

where M P \ ~ 2.4 • 10 18 GeV is reduced Planck scale, which will be treated as natural cut off for 
all higher dimensional non-renormalizable operators. In our approach, top quark (and possibly 
bottom quark and tau lepton, in case of large tan f$) will get mass at renormalizable level. Yukawa 
couplings of light families emerge after U(l)p flavor symmetry breaking. Thus, the hierarchies 
between Yukawa couplings and CKM mixing angles will be expressed by powers of small parameters 
e,e< 1. 

Yukawa textures 

In the charge assignment we need to fix the values of a, b and c. If their ratios remain arbitrary 
there will be more than one extra U(l) symmetry, and that we have to avoid. Together with 
fixing a,b,c, the values of p,q and k (in (5) and (6)) should be selected in such a way as to have 
phenomeno logically viable quark and lepton Yukawa textures. It turns out, that one selection 
leading to attractive Yukawa sector, is the following: 

{«,&,c} = {-^^|} > P = q = k = 0. (9) 

In this case, in Eq. (6) we pick up / = 3, which fixes charges of H,H as (Qh,QhY =3 — (0,2). 
The charges of matter (10 + 5)-plets (and also SU(5) singlets) are also fixed and we will make 
the following identification: Q Wi = {2,-1,0} , Qi % = {0,-3,-2}, where % = 1,2,3 labels the 
flavor. The model's field content and corresponding U(Y)f charges are given in Table 1. With this 
assignment, 10 • 10if and 10 • 5^-type Yukawa couplings are given by 




10! 10 2 
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e 2 e \H 
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5 2 


5 3 


10i 
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io 2 




e 2 




io 3 




e 





(10) 



3 Details of symmetry breaking, with general values of (X) and (X), is given in Ref. [8]. 



4 



Table 1: U(1)f charge assignment for the model's states. 
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where in front of each entry, dimensionless couplings(~ 1/5 — 5) are assumed. As we will see shortly, 
good fit is achieved for e ~ 1/10 , e ~ (0.05 — 0.2)e 2 . On the other hand, with these values, the 
matrix elements (1,1), (1, 3), (3, 1) are so suppressed, that they are irrelevant and we can set them 
equal to zero. Thus, for all practical purposes, we can investigate the Yukawa matrices: 



Yu,d,e oc 



with zero textures. 



(° 

e 



e 

^2 



(11) 



Quark masses and mixings 

Using the basis q T Yuu c h u and q T Y D d c hd, without loss of generality we can parameterize up and 
down Yukawa matrices at GUT scale to have forms: 
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(1 ^p^M 





Ae 2 
a d e 2 e^ d 
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(12) 



(13) 



We have made field phase redefinitions in such a way that, in this basis, CKM matrix remains 
unity and in Yu only one phase £ M appears. The phases ip and ip' will not contribute to the quark 
masses, but will be important for the CKM matrix elements. From Eqs. (12) and (13), in a fairly 
good approximation we obtain the following relations valid at GUT scale: 



c ~e 2 k,e 2 e* 



As _ 2 \a d e 2 e^ d -bb'\ 

y h - e i z (uz\2 



cA — \ — - 



Ae 2 ^ 




(1 + (b'e) 2 ) 



2\3/4 



(14) 



(15) 



v fe i + {vzy 

These relations help to find a good fit. With proper selection of input parameters e, A, a Ut d, 
b, b', c, c', k, £ U! d, (p, p' we can get desirable values for fermion mass hierarchies and CKM mixing 
angles at GUT scale. Then, using RG we can calculate these ratios at low scales: 
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A, 



3 ^ u 

VtVb 
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V a p\ Mz = VtVb V a p\ Ma , if (a/3) = (ub,cb,td,ts) 

V a p\ Mz = V <*p\ Mg > if ( a P) = (ud,us,cd,cs,tb) , (16) 

where RG factors 

* = ex p^y mt > * = ^(jeW™, A ^ ln/i J ' * = exp v^l, A ^ ln/i J (17) 

are given in 1-loop approximation. 

We will consider two cases with low/moderate and large values of the MSSM parameter tan/3. 

Fit for tan (5 = 5 - 15 

We take experimental value m t (m t ) = 163.68 GeV, determining top Yukawa coupling at weak scale, 
and with tan/3 = 5 — 15 we find rj t = 1.097, r\ b ~ r\ T ~ 1 . For this case, good fit is obtained for the 
following values of input parameters: 

e = 0.0847 , A = 0.476 , a u = 0.6 , a d = 3.7 , 

b = _0.798, 6' = -7.14, c = 0.037, fc = 0.864 , 

£ u = o , £ d = -0.065 , up = -2.696 , <p' = -0.97 . (18) 
These at GUT scale give 

at /i = Mq : ^ = 5.51 • KT 6 , ^ = 0.002835 , ^ = 5.64 • 10~ 4 , ^ = 0.0111 , 
At At \ b X b 

1^1 = 0.2243, 1^1= 0.0383, |K&| = 0.00318 , p = 0.118, = 0.34 , 

where p + irj = - \ U r d Xh ■ 

Performing renormalization (using (16) and [13]), at low scales we get (with input m t {m t ) = 
163.68 GeV, m b (m b ) = 4.24 GeV): 

(m u , m d , m s , m c ) (2 GeV) = (2.1, 4.64, 91.69, 1082) MeV 

at p = M z : \V US \ = 0.2243 , \V cb \ = 0.042 , \V ub \ = 0.00349 , p = 0.118 , rj = 0.34 . (19) 
These values of masses and CKM matrix elements are in good agreement with experiments [14], [15]. 

Fit for tan f3 = 55 

In this case we have r\ t = 1.114, rj b = 1.158, i] T = 1.105 . Input parameters are selected as: 

e = 0.0723 , A = 0.53 , a u = 0.545 , a d = 6.81 , 

b= -0.777, b' = -11.58, c = 0.0375, k = 0.783 , 
£ u = -0.055 , U = -0.0593 , ip = -2.73 , up = -0.98 , (20) 
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giving at GUT scale 

at n = M G : ^ = 4.5 • 1(T 6 , ^ = 0.002367 , ^ = 3.73 • 10~ 4 , ^ = 0.00723 , 

At \t A5 At 

\V US \ = 0.2259 , \V cb \ = 0.0317 , \V ub \ = 0.00272 , p = 0.135 , rj = 0.345 . 

The renormalization procedure (using (16) and [13]) gives at low scales (with input m t {m t ) = 
163.68 GeV, m b {m b ) = 4.24 GeV): 

(m u , m d , m„ m c ) (2 GeV) = (2.08, 4.85, 93.86, 1096) MeV 

at /i = M z : \V US \ = 0.2259 , \V cb \ = 0.0409 , \V ub \ = 0.00351 , p = 0.135 , rj = 0.345 . (21) 

These agree well with experiments. 

Charged lepton sector 

Now let us discuss the charged lepton sector. Relevant Yukawa couplings originate from 10 • 5 ■ H- 
type interactions of Eq. (10) (while in practice Y E has the structure of Eq. (11)). Without breaking 
the SU(5) symmetry in these interactions, one would get the asymptotic relation M D = M|, which 
is unacceptable and is a well known problem for minimal SUSY SU(5) GUT. However, by some 
specific extension, care can be exercised to solve this problem [16]. Without specifying origin of 
SU(5) breaking in this sector, we assume that it happens (i.e. SU(5) symmetry breaking) in the 
sector of light families. Thus, in analogy of Y D (see Eq. (13)), in a basis l T Y E e c h d) we parameterize 
Y E to have the following form 

/ k e X e e 2 \ 

Y E ~ \ e e 2 k 22 a d e 2 e^ b'e A° , with A° = A° . (22) 
\ be I J 

In Y E only one complex phase £ e appears. Remaining phases are rotated away by proper phase 
redefinitions of the / and e c states. With {/C22, k e , A e , ^ e } 7^ {l,k, A,^} we can avoid the relation 
rrT 1 = m^' wnne keeping m" = mj (at the GUT scale). Good fit can be obtained with 

for tan/3 = 5-15, A e = 2.51 , k e = 0.082 , k 22 = 4.517 , £ e = -0.065, 

for tan/3 = 55, A e = 3.11 , k e = 0.07656 , k 22 = 3.3677 , £ e = -0.06 , (23) 

and remaining parameters given in Eqs. (18) and (20) respectively. With these we obtain 

at /i = Mq , for tan/3 = 5- 15, ^ = 2.787 • 10~ 4 , ^ = 0.05883, 

at n = M G , for tan/3 = 55, ^ = 2.065 • 10" 4 , ^ = 0.0436. (24) 

These lead to 

m e (m e ) = 0.511 MeV, m M (m M ) = 105.66 MeV , m T (m T ) = 1.777 GeV , (25) 

in agreement with experiments. The mixing angles originating from the charged lepton sector, for 
tan/3 = 5-15, are {0| 3 , 6 e l2 , 6 e n } ~ {31.5°, 1.02°, 0.61°}. While for tan/3 = 55 we got {0f 3 , 0f 2 , ^ 3 } ~ 
{40.2°, 0.93°, 0.78°}. Note that while 0| 3 is large (but not sufficiently), the and 6 e n are too small. 
This means that neutrino sector should be responsible for generating proper values of the lepton 
mixing angles. 
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4 Neutrino Sector 



To build the realistic neutrino sector, we apply the singlet states 1.1,2,3 (with U(1)f charges given in 
Table 1) as right handed neutrinos. Their Dirac type couplings (to 5j states) and the mass matrix 
respectively are given by: 
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5 3 
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e 2 
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e 2 



(26) 



where M* is some mass scale and in the entries of these matrices the dimensionless couplings are 
omitted. Integration of heavy 1, states leads to 3 x 3 mass matrix for the light neutrinos: 



(e 2 ee 2 ee 
ee 2 a 2 e 4 a(3e 3 ) m , ^ 21 > 
ee afc 3 (3 2 e 2 

with fh ~ ^ 2 and a, /3 are some dimensionless couplings. Note that, M^'s 2 — 3 block's determi- 
nant is zero. It is convenient to work in a basis where charged lepton mass matrix is diagonal, i.e. 
rotate whole lepton doublets by unitary matrix which diagonalizes the matrix Y E Y^. In this basis, 
the weak leptonic current is diagonal and the neutrino mass matrix can be denoted by M v . The 
convenience of this basis is that the diagonalizing matrix U : 

U T M V U = Mf 8 (28) 

will coincide with the lepton mixing matrix. The latter, in a standard parametrization, has the 
form: 

(C13C12 C13S12 s 13 e- tS \ 

-c 23 si 2 - s 23 Si3Ci 2 e i<5 c 23 ci2 - s 2 3Si3Si 2 e l<5 S23C13 P2 (29) 
523^12 - C 23 Sl3Cl2e l5 -S23C12 - C 2 3Si3S 12 e lS C 23 C 13 ) 

with Sij = sin 9ij and c^- = cos%. The phase matrices Pi ;2 are given by: 

Pi = Diag (e" iwi , e~ iul2 , e"^ a ) , P 2 = Diag (l , e~ ipl/2 , e~ ip2/2 ) , (30) 

where ^1,2,3, Pi, 2 are some phases. 

To get some feeling about the results, obtained from the neutrino mass matrix, let us first ignore 
9f 2 and 9f 3 mixings. Since these angles are small, the picture qualitatively will remain unchanged. 
(Effects of these mixing angles are discussed in detail in an Appendix). With Q\ 2 -,Q\ 3 1, in a 
good approximation M v can be written as: 
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b 2 


ab 
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ab 


a 2 



(31) 
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Figure 1: Region (i): Needed values of K, realizing normal hierarchical neutrino masses. Region 
(ii): Values of K within considered scenario with normal ordering of neutrino masses of Eq. (56). 



Note that 2 — 3 block's determinant of the matrix (31) is also zero: 
This, using (29), leads to the following interesting relation 4 



^(2,2)^(3,3) _ ( ^(2,3) )2 



0. 



tan 2 6» 



13 



m 3 
m 2 



<? 2 e 



mi 



.12 



(32) 



Using recent results from the neutrino experiments [1-3], we can easily verify that the relation of 
Eq. (32) is incompatible with normal hierarchical neutrino masses. This conclusion remains robust 
taking into account the effects of 1 — 2 and 1 — 3 rotations coming from the charged lepton sector. 
With these, instead of Eq. (32) we have the exact expression 



tan 2 6» 



13 



7^3 

m 2 



\s 12 e 



+ ^c 2 I 
+ m 1 Cl2 ' 



+ 



K 2 



17111712 



(33) 



where K is real and k is some phase. Derivation of (33) and forms of K, k are given in Appendix 
(see Eqs. (53)- (55)). One can investigate for what values of K, desirable values of 9 13 are obtained. 
In Fig. 1, region (i) corresponds to the values of K as a function of 777,3, which give #13 ~ 8.9°. On 
the other hand, region (ii) shows values of K obtained within considered scenario (for 9\ 2 = 0.016 
and 8f 3 = 0.0136). We see that points of region (ii) are well below from points of region (i). While 
Fig. 1 corresponds to the best fit values of the neutrino oscillation parameters [1], the conclusion 
is same by taking them within 8a error bars. This demonstrates that within considered model, the 
normal hierarchical neutrino mass scenario can not be realized. 

On the other hand, inverted hierarchy in neutrino masses is possible within considered SU (5) x 
U{1)f model. 5 This is demonstrated in Fig. 2. Green dashed region includes points captured 



4 See Appendix for exact expression, detailed derivation and related discussion. 

5 Worth pointing that it is usually hard to get inverted neutrino mass scenario within GUTs [17]. See however [18] 
(within SO(10)) and [19] (within Eq GUT) with inverted hierarchical neutrino masses. 
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0.001 0.002 0.003 0.004 0.005 



m 3 (eV) 

Figure 2: Inverted hierarchical neutrino mass scenario. Green dashed region shows allowed values 
of (7713, #13). Red bold curves (a) and (b) represent the dependance of #13 on m 3 for phases 
{pi, 6 + P2, UJ3 — 0J2, K } — {0, 0, 0, 0} and {tc, tc, 0, tc} respectively. Curve (a) gives largest possible 
values of #13, while curve (b) - lowest ones. Two horizontal lines are upper and low experimental 
bounds of 6*13 within the la. Dashed curves would had been obtained (instead of bold ones) with 
K -»■ in Eq. (33). 



by two border bold curves (obtained via Eq. (33) within our model) and two horizontal lines 
(corresponding to the experimental values of #13 within la). This figure corresponds to the best 
fit values of #12, #23, ^ m 1 \ an d Am^ tm , while free phases (see Eqs. (54), (55)) are varied within 
full ranges. Dashed lines correspond to the case with K — > 0. Thus, inclusion of the charged 
lepton sector somewhat extents the allowed region. All this demonstrates that inverted hierarchical 
scenario is easily realized. Fig. 2 shows that, the allowed region for m 3 is fixed as: 

0.0008 eV < m 3 < 0.0044 eV , (34) 

and using (58) we get: 

1/2 



m 1 ~ 0.04852 eV x ^1 + (- 
m 2 ~ 0.0493 eV x ^1 + (— 



m 3 



04852 eV 



. , 1/2 
m 3 - x 



0493 eV 



(35) 



These imply Yl m i ~ 0-1 satisfying the current bound [20] obtained from cosmology. Moreover, 
for neutrino less double /3-decay parameter mpp = \ ^ U^m^ we obtain: 

% ^ \(?\2 m \ + s 2 12 m 2 e' lpl I , (36) 

leading to: 

0.011 eV < m w < 0.05 eV . (37) 
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Future experiments will be able to test viability of this scenario [21]. 

In summary, we have presented supersymmetric SU(5) GUT supplemented with non-anomalous 
U(l)p flavor symmetry. Anomaly cancellation condition restricted the field content (dictated three 
families of 10 + 5 matter), as well as U(1) F charge assignment. Texture zero Yukawa matrices 
gave natural understanding of hierarchies between charged fermion mass and mixings. Model 
automatically involves SU (5) singlet states utilized as right handed neutrinos. Inverted hierarchical 
neutrino mass scenario is predicted within considered model. Other phenomenological issues, such 
as doublet-triplet splitting, proton decay etc., left beyond the scope of this paper, will be addressed 
elsewhere within more general class of models [12] supplemented by anomaly free U(1)f symmetry 
[8]. 
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Appendix: Effects of 0?- on the Neutrino Sector 

In this appendix we work out the details of contributions from the charged lepton sector to the 
neutrino sector. In particular, as was pointed out in Sect. 4, we study impact of 6\ 2 an d #13 mixing 
angles. With this study, we prove that within considered SUSY SU(5) x U(1) F scenario, only 
inverted hierarchical neutrino mass scenario is realized. 

Charged lepton mass terms e T M E e c get diagonalized by transformations e = L* e e' , e c = R e e c/ , 
where L e , R e are unitary matrices such that 

L\M E R e = Mf s . (38) 

Let us rotate the neutrino states v by the same unitary transformation as e-states: v = L* e v' . With 
this, the weak current remains diagonal: ej^u = e'j^v' . On the other hand, the neutrino mass 
couplings \v T M v v become \v' T M v v' with 

M v = L\M u L* e ■ (39) 

Upon transformation v' — U v" , the neutrino couplings can be diagonalized, i.e. 

U T M U U = M„ Diag , (40) 

and finally the weak current will be e'^^Uu" . Thus, the matrix U in (40) coincides with the lepton 
mixing matrix. 

From Eqs. (39) and (40) we obtain 

L\M v L* e = U*M^ iag U^ . (41) 
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Unitary matrix L e can be written as 

L e = P l x L 23 L l3 L 12 P l 2 , (42) 

where P[ 2 are some diagonal phase matrices, L± 2 , L 23 and L i3 correspond to the rotation angles 
9f 2 , 9 23 and 9{ 3 respectively. Without loss of generality, the unitary matrices L 12 and L 23 can be 
taken to be real orthogonal matrices. Since within our scenario 9\ 2 and 9\ 3 are small, we can write 



L 13 L 12 ~ P/(l + V)Pl' , (43) 



where V is real: 





f ° e 


6 12 


6 13 














\ A 13 









, with s{ 2 = sin 0^ , s^ 3 = sin 0f 3 . (44) 

Without restricting any generality, we can take P\ = Pi* and using (42), (43) in (41), we obtain 

(1 + T T )M* u (l + T) = UM™* S U T , (45) 

with 

M v = P;*Ll 3 Pl*M u P*L* 23 P' L * . (46) 

Because of smallness of s\ 2 and sf 3 , further we use approximation and keep first powers of these 
angles (and thus first powers of the matrix T). With this, from (45) we get 

Ml + Y T M* V + M;T = T , with T = UM^ iag U T . (47) 

Note, that since 2 — 3 block's determinant of M v is zero, similar applies to the 2 — 3 block of the 
matrix M v (see Eq. (46)). Thus, M* can be parameterized as 



e 


c 


d 


c 


b 2 


db 


d 


db 


a 2 



M*=\c b 2 ab I . (48 
\ d db 

With this, using matrix relation in Eq. (47), we derive 

e - 2(cs e 12 + ds e 13 ) = T u , c + es e 12 - b(bs e 12 + as e 13 ) = T 12 = T 21 , 

b 2 + 2cs{ 2 = T 22 , d + e Sl e 2 - d{bs{ 2 + a Sl e 3 ) = T 13 = T 31 , 

a 2 + 2ds e 13 = T 33 , a& + cs e 13 + rf Sl e 2 = T 23 = T 32 . (49) 

By iteration (keeping 0(s{ 2 ) and 0{s{ 3 )) we obtain from (49): 

d 2 = T 33 — 2s^ 3 T 13 , b = T 22 — 2s\ 2 T\ 2 , 

a6 = T 23 - s e 12 T 13 - Sl e 3 T 12 . (50) 
These three expressions, by eliminating d and b, give the following relation: 

T 23 — T 22 T 33 = 2s^ 2 (T 23 T 13 — T 12 T 33 ) + 2s^ 3 (T 23 T 12 — T 22 T 13 ) . (51) 
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Substituting Tij elements (see Eq. (47)) in (51) we obtain 
m 1 m 2 (U2iU 3 2-U22U 3 i) 2 + m 1 m 3 {U2iU 33 -U2 3 U 3l ) 2 + m2m 3 {U22U 33 -U2 3 U 3 2) 2 =-2s e l2 K l -2s e lz K2 (52) 
with 

K x = \ ( mi m2 sm26 13 S2 3 e t{5+p2) + (m ie ipl -m 2 )m 3 sin2^ 12 c 23 ) e -^+P^+^+^) , 

K 2 =- {m x m 2 sm26 l3 C2 3 e i{&+p2) - [m x e ipx -m 2 )m 3 sin2^ 12 s 23 ) e - i ^ 1+ " 2+wl + 2w2 + W3 ) . (53) 

Z 

Using the form of C/ of Eq. (29) in left hand side of (52), after some simplifications we obtain 

- tan 2 e 13 e i{ - 25+ ri = s\ 2 e^+ — c 2 2 V (s^K.+ s^K,) e i ^+f»+*»+*>*) . (54) 

m 2 \ m 1 J mim 2 cf 3 

Introducing notations 

K 2 = -J- Ist^+s^l , 

c 13 

k = Pi + p 2 + 2w 2 + 2w 2 + Arg (s^Kx+st^) - Arg (mis? 2 e <pi + m 2 c 2 2 ) , (55) 

from (54) we get Eq. (33) - the expression for tan 2 6 13 . 

Having (54), we can now examine possibilities of realizing normal and inverted hierarchical 
neutrino mass scenarios within our model. 

Excluding normal hierarchical neutrino mass scenario 

Let us first see if normal hierarchical neutrino masses are possible. In this case, m 3 > m 2 > mi and 
observed mass squire differences are Am 2 ol = m\ — m\ and Am 2 tm = m\ — m\. Thus, two masses, 
say mi and m 2 , can be expressed as 

mi = \jm\- Am s 2 ol - Am 2 tm , m 2 = \jm\- Am 2 tm . (56) 

While Am 2 ol = m 2 — m 2 , Am 2 tm = m\ — m 2 are both measured, the m 3 is unknown yet and we will 
treat it as a free parameter. Taking into account (52), we can easily verify that m x = is excluded. 
Thus, m 3 > a/ Am 2 ol + Am 2 tm ~ 0.05 eV. On the other hand, we can also have an upper bound 
for m 3 , set from the cosmological bound on a sum of three neutrino masses J^mj < 1 eV [20,21]. 
This, taking into account (56), gives m 3 < 0.34 eV. Therefore, we will vary m 3 in a range 

0.05 eV < m 3 < 0.34 eV . (57) 

With help of (33) and using the best fit values of quantities Am 2 ol , Am 2 tm and neutrino mixing 
angles [1], we can see what values of K are needed. In Fig. 1, dashed region (i) represents such 
values of K (versus m 3 ). For fixed value of m 3 , the multiple values of K are obtained because 
of free phases appearing in (33). On the other hand, with (55) and (53) within our model with 
sf 2 = 0.016, sf 3 = 0.0136 we can calculate K for different m 3 and remaining phases. Region 
(ii) of Fig. 1 corresponds to this. We see that regions (i) and (ii) do not overlap and therefore 
conclude that normal hierarchical neutrino mass scenario is not realized within considered model. 
This conclusion remains robust even varying the values of Am 2 ol , Am 2 tm , #i 2 , # 23 , #i 3 within 8cr 
error bars. 
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Compatibility with inverted hierarchical neutrino masses 



As turns out, the inverted hierarchical neutrino masses blend well with relation (33). In this case, 
Amg 0l = m\ — m\ and Am^ tm = m\ — m\ and thus: 



If we set K — > in (33), we can easily see that for certain values of m 3 and p± all observable can be 
obtained within experimentally preferred ranges. Inclusion of K do not change this positive result, 
but just offers slightly different choices of m% and various phases. For illustration see Fig. 2, with 
corresponding discussion starting in a paragraph right before Eq. (34). 
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